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Abstract 

Q This article explains why a paper by Heinz G. Helfenstein entitled "Ovals with 

ch equichordal points", published in J. London Math. Soc. 31, 54-57, 1956, is incorrect. 

We point out a computational error which renders his conclusions invalid. More im- 
portantly, we explain that the method cannot be used to solve the equichordal point 
problem with the method presented there. Today, there is a solution to the problem: 
Marek R. Rychlik, "A complete solution to the equichordal point problem of Fujiwara, 
Blaschke, Rothe and Weizenbock", Inventiones Mathematicae 129 (1), 141-212, 1997. 
However, some mathematicians still point to Helfenstein's paper as a plausible path 
to a simpler solution. We show that Helfenstein's method cannot be salvaged. The 
fact that Helfenstein's argument is not correct was known to Wirsing, but he did not 
explicitly point out the error. This article points out the error and the reasons for the 
failure of Helfenstein's approach in an accessible, and hopefully enjoyable way. 

CM 

• • 

1 The Equichordal Point Problem 

The equichordal point problem enjoyed significant popularity since its original formulation by 
Fujiwara in 1916 and Blaschke, Rothe and Weizenbock in 1917, because it can be formulated 
in easy to understand terms of elementary geometry, and it is hard to solve. The starting 
point is the definition of an equichordal point: 

Definition 1. Let C be a Jordan curve and let O be a point inside it. This point is called 
equichordal if every chord of C through this point has the same length. 

Then the equichordal point problem may be formulated as follows: 

Question 1. Is there a curve with two equichordal points? 
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Why two? Because circles and a lot of other shapes have one equichordal point, and 
because Fujiwara pointed out that it is impossible for a shape to have three equichordal 
points. 

The full solution to the problem appears in the article [lj. The paper is considered (and 
it is!) rather hard to read and its length is 72 pages. Thus, although it is a great resource for 
anyone studying this and related problems, it is not always easy to extract the information 
one needs. 

In this set of notes we use the information provided in [1] to construct a counterexample 
to a published article by Helfenstein [2] in 1956. Helfenstein made a claim which would 
lead to a simple solution of the equichordal point problem (under 10 pages, perhaps) if it is 
augmented with a few relatively easy facts to prove. It has been hope of many that such a 
simple proof exists. However, as we will see, Helfenstein's paper is incorrect, and thus there 
is no hope for a simple proof, at least along the lines of Helfenstein's argument. 

The convex geometry community, in which the equichordal point problem, and our solu- 
tion of it, are quite well known. The community has had an especially hard time coming to 
grips with the hard analytical methods used in our article (and also prior articles of Wirs- 
ing [7] and Shafke and Volkmer [3]). We found on several occasions that the argument of 
Helfenstein continues to have some legitimacy because no one has explicitly shown where it 
is incorrect [5]. At the end of this paper we cite Grunbaums's argument from |5] which is 
representative of this opinion, although the experts on the problem (including Wirsing) have 
clearly dismissed Helfenstein's paper. Therefore, it will be beneficial to analyze Helfenstein's 
argument from today's perspective, and explain why it is incorrect. We hope that the read- 
ing is entertaining and allows one to understand some of the trappings of the problem, and 
perhaps even appreciate the length and complexity of our solution. 

Helfenstein's paper contains an incorrect statement which must have resulted from an 
error in a mundane calculation, involving only elementary calculus. This will be clear from 
what follows. With the aid of a Computer Algebra System (CAS), we reconstructed and 
corrected the intermediate calculations, and arrived at the opposite conclusion, which clearly 
shows an error in Helfenstein's argument in an elementary way. 

More importantly, the main idea of Helfenstein's paper is also incorrect, and it cannot 
be salvaged by simply correcting the error in calculation he made. We show this in the 
strongest possible way: we construct a counterexample by referring to the relevant portions 
of [lj. However, we will make the argument as simple and as self-contained as possible. 

2 A summary of Helfenstein's paper 

In 1954 Heinz Helfenstein submitted an article [2], in which he claims that there is no oval 
with 2 equichordal points that is 6 times differentiable. He calls a curve with 2 equichordal 
points a 2e-curve. We will use this abbreviated term below, as synonymous with "curve with 
2 equichordal points". 

We proceed to summarize the terminology, technique and results of Helfenstein's paper. 
For simplicity, we will assume that the curve C under consideration is a convex oval and it 
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Figure 1: Illustration of Helfenstein's notation. 

is symmetric in various ways (cf. Figure [T]): 

1. It is symmetric with respect to a point O inside C . 

2. Let R and S be the two hypothetical equichordal points. We may assume that O 
bisects the interval RS. 

3. C is symmetric with respect to the straight line passing through R and S. 

4. Let A and B be the two points collinear with R and S which belong to C. We may 
assume that the distance \AB\ = 1. Moreover, due to prior symmetries, O bisects AB. 

5. The distance BS is called c. By the symmetries assumed, we have AR = c. 

6. The construction is valid when < c < 1. The order of the points A, R, O, S and B 
on the straight line on which all these points lie is: 

(a) A,R,0,S,B if < c< 1/2; 

(b) A, S, O, R, B if 1/2 < c < 1. 

For c = 1/2, R = = S. 

It should be stated that the above picture of a hypothetical 2e-curve is correct, based on 
many independent analyses. In particular, the symmetries are well established. 

The next assumption used in the paper [2] is that locally near A and B the curve C 
may be represented by a graph of a function. Helfenstein uses two orthogonal coordinate 
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systems, one centered at A and one centered at B. The coordinates of the system centered 
at A are called £ and rj, and the positive direction of the r] axis is AB. The coordinates of 
the system centered at B are called x and y, and the positive direction of the y axis is BA, 
so that the y and rj axes point in the opposite directions. We assume that C near the points 
B and A may be represented by the equations: y = f(x) and rj = /(£), respectively. There 
is an agreement of results supporting the claim that C is represented by a smooth function 
f(x) near the points A and B. It should be emphasized that f(x) is locally defined, i.e. its 
domain is some interval (— e, e), where e > 0. There is a proof that e may be as large as 1/2, 
but this will not be material in these notes. One can also prove that f(x) is real- analytic, 
i.e. it may be represented by a power series convergent on the interval \x\ < e. Again, the 
analyticity is not material in these notes, but Helfenstein assumes that the function has six 
derivatives. It should be noted that Helfenstein's assumption in regard to differentiablility 
is faulty. The fact that he assumes six-fold differentiability is a result of a computational 
error, as will be demonstrated below. 

The next important construction in Helfenstein's paper is that of a functional equation 
satisfied by f(x). The derivation presented in the Helfenstein's paper is correct, and is 
consistent with the construction used in our solution of the equichordal point problem pp. 
It nicely illustrates the transition from geometry to analysis, which is a hallmark feature of 
the equichordal point problem. We will repeat the Helfenstein's construction here. 

Let P(£, 77) be a point near A and let Q(x,y) be a point near B, both on the curve C 
and both represented in the respective coordinate systems. Let P' and Q' be the orthogonal 
projections of P and Q onto the line AB. Let T be the projection of P onto the line QQ', 
perpendicular to AB. Helfenstein observes that the triangles QQ' S and QTP are similar. 
From this observation, the following equations result: 

x x + £ 



\J x 2 + (c — y) 2 1 

c-y l-y-rj 



a/x 2 + (c-y) 2 

By solving with respect to £ and 77 we obtain: 



c X 
£ = , — x 

a/x 2 + (c-y) 2 



y' 2 

( c -y) 



+ (!-</)■ (2) 



a/x 2 + (c-y) 2 

By construction, y = f(x) and rj = /(£). Therefore, we obtain this functional equation: 

f ( x - x\ = ~ (c ~ /(x)) _ + a - fix)) (3) 

\y/x* + {c-f{z))* J ^x 2 + (c-f(x)) 2 

The manner in which Helfenstein uses this equation is also well established: we repeatedly 
differentiate both sides at x = to obtain the consecutive derivatives of f{x) at x = 0. Thus, 
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we try to solve the equation by finding the /(x)'s Taylor series at x = 0. Helfenstein uses the 
following notation, defining coefficients of the Taylor expansion at x = 0: for n — 0, 1, 2, . . .: 



d n f(x) 



n\ a 



(4) 



x=0 



Often we will write a n (c) instead of a n when it is necessary to consider the dependence of 
a n upon the parameter c. 

The consensus of several methods is that these equations can be used to determine a n 
by a recurrence relation, and thus determine the Taylor expansion of f(x) at x = up to an 
arbitrary order. Moreover, the symmetries imply that f(x) is an even function: 



This implies that a n = for odd n. Moreover, /(0) = follows from the assumptions made, 
that C passes through A and B. 

We come to a point where Helfenstein makes a calculation error in calculating the third 
non-trivial coefficient, clq. Only simple calculus (chain rule) is involved. Calculating ci2 and 
a 4 by hand would test anyone's patience, but today is conveniently done with the aid of a 
Computer Algebra System (CAS). Helfenstein calculated 02 and 0,4 correctly. Calculation of 
a§ must have been very challenging without a CAS, and indeed it resulted in an important 
error which affects the entire argument. 

We wrote a simple CAS program which determines the coefficients a n for n up to 10. In 
theory, the program can find a n for arbitrarily large n, but even CAS consumes an amount 
of time that probably grows exponentially with n. The CAS we used is the open source, 
free system Maxima [6j, although any other CAS can solve this problem. We included our 
program as Appendix A. 

The results are presented below for even n only. Moreover, we list numbers 



which are analytic in z iff a n are invariant under the substitution c4 1 — c. Thus, it is much 
easier to read off the invariance by looking at h n . 

The program generated the coefficients in standard TgX format. Additionally, the pro- 
gram factored a n as rational functions of c, for easy comparison of 02 and 04 with Helfenstein's 
paper. Here is the result: 



2 (2 c 2 - 2 c + 1) 

12c 4 - 24c 3 + 12c 2 -1 

8 (2c 2 - 2c +1) 2 (2c 4 - 4c 3 + 6c 2 - 4c +1) 
80 c 10 - 400 c 9 + 680 c 8 - 320 c 7 - 500 c 6 + 940 c 5 - 712 c 4 + 284 c 3 - 52 c 2 + 1 

16 (2c 2 -2c+l) 4 (c 4 -2c 3 + 5c 2 -4c+l) (2 c 4 - 4 c 3 + 6 c 2 - 4c + 1) 



f(-x) = f{x). 




(5) 



1 



5 



and 



b 2 

h 



z + l 



3 z 2 - 6 z - 1 



z* + 8^ 3 + 14^ 2 + 8,2 + 1 

10 z 5 - 110 z 4 - 20 



204 z 2 + 42z + 2 



z 8 + 24 + 172 z 6 + 488 z 5 + 678 ^ 4 + 488 z 3 + 172 ^ 2 + 24 z + 1 



On the web we listed the coefficients up to a w without folding or breaking them up. The 
corresponding b n are clearly analytic, i.e. do not contain half-integer powers of z. This 
means that a n is invariant under the substitution c >->■ 1 — c for n up to 10. We can push 
this calculation further, up to, say, n = 20, always with the same result: it is invariant under 
this substitution. 

Helfenstein's paper contains correct expressions for a 2 and a^. However, he did not 
include the expression for a 6 . Since he derived a false conclusion about a 6 , as we will see 
below the only possible explanation is that he made a computational error in the intermediate 
calculations. Helfenstein's argument is founded on an unproven claim that if a 2e-curve exists 
for some value c then a n must be invariant under the substitution 1 — c. More precisely, 
if we consider a n = a n (c) (i.e. as a function of c) then the condition a n (c) = a n (l — c) is 
necessary (according to Helfenstein) for a 2e-curve to exist for a particular value of c. 

Finally, we are ready to explain Helfenstein's main argument, and how he arrived at the 
erroneous six-fold differentiability condition. He correctly noted that the expressions for a 2 
and CZ4 are invariant under the substitution c i— > 1 — c. He then considers ciq as a candidate 
for a coefficient which is not invariant under this substitution. In contradiction with our 
findings, he claims that it is not invariant under the substitution c 4 1 — c. Helfenstein 
writes: 

A sixth differentiation finally yields an expression for a 6 (c) which is not identical 
to a 6 (l — c). 

The form of a§ is omitted in Helfenstein's paper and the intermediate calculations of it are 
missing. He then proceeds to determine specific values of c which solve the equation: 

a 6 (c) = a 6 (l - c) 

He claims that the above equation is equivalent to a certain polynomial equation of the 9-th 
degree: 

144 c 9 - 648 c 8 + 1176 c 7 - 1092 c 6 + 168 c 5 + 798 c 4 - 846 c 3 + 357 c 2 - 59 c + 1 = 0. 
Subsequently, he demonstrates that equation does not have roots in the range 

2-y/3 2 + V3 
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which is known to be the region of c, outside of which there is no 2e-curve, based on ele- 
mentary arguments which preceded Helfenstein's paper. Clearly, the 9-th degree polynomial 
and the subsequent conclusions are a result of a calculation error. 

Helfenstein's claim is that there are no 2e-curves which are six-fold differentiable. The 
reason is that he needs this much differentiability to calculate a 6 . As we have shown, the six- 
fold differentiability of the function f(x) at x = is not sufficient to disprove the existence 
of a 2e-curve. Moreover, we verified with a CAS that Helfenstein's argument fails for curves 
which are ten-fold differentiable, on the basis of our calculations of a n and b n for n up to 10. 

Of course, now it is time to pause, and suggest that the argument fails for all n. This 
shows that Helfenstein's approach cannot succeed, even if we had unlimited computing power 
and find as many coefficients a n as necessary. 

3 More differentiability does not help 

One could hope that by finding more derivatives of f(x) we will eventually find a coefficient 
a n (c) for which a n (c) and a n (l — c) do not coincide. In this section we will show that a n (c) = 
a n (l — c) for all n, given that all derivatives up to n exist. This, of course, demonstrates 
that Helfenstein's method cannot disprove the existence of a 2e-curve. 

The key point is to understand the local existence and uniqueness of solutions of Helfen- 
stein's functional equation. It should be noted that his paper is an attempt to disprove 
local existence. What he missed is the fact that the locally defined solution to his functional 
equation does exist! Furthermore, he missed that the local existence does not imply that a 
2-e curve exists. 

The existence and uniqueness results are contained in the Inventiones article pQ, but we 
will formulate those results here in a manner more suitable for these notes. 

The graph of a solution to the Helfenstein's functional equation gives rise to two Jordan 
arcs, Ca and Cb, contained in a neighborhoods of A and B respectively. The arcs Ca and 
Cb are defined by the equations 

v = /(£), 
y = f(x) 

in the respective coordinate system. A picture is worth a thousand words. Thus, if the 
reader still cannot imagine how the two arcs Ca and Cb may connect together when they 
are maximally extended, without forming a 2e-curve, a plausible scenario can be visualized 
by a schematic drawing in Figure [2j The alternative would be for the two arcs to meet 
exactly and form one smooth curve. The fact that they do not meet in this manner is the 
thrust of the Inventiones article [IJ . It should be noted that the figure has perfect reflectional 
symmetry, both with respect to reflections in the line AB and in the center point O. Because 
Helfenstein did not look away from the line AB (i.e. outside the "cuts", which mark the 
ends of Jordan arcs Ca and Cb) he failed to notice that his assumptions may be satisfied by 
a local equichordal configuration. And indeed, this is what really happens. While the above 
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Figure 2: A schematic figure of local curves connecting. 



figure is only a schematic, we and others performed numerical computations which are in 
perfect agreement with the above picture in regard to its topology. It should be noted that 
the oscillations by the curves near A and B continue up to the line AB. The size of the 
oscillations settles down to a fixed, positive amplitude. 

The extremely important idea in understanding the equichordal point problem has been 
that the problem should be phrased as a problem bout iterations of mappings, i.e. should 
be framed as a problem of dynamical systems theory. To remain faithful to Helfenstein's 
notation, we define a map G c : U c — > R 2 , where U c C R 3 is an open, punctured unit disk 
centered at S(0, c): 

U c = {(x,y) Gl 2 : < x 2 + (y - cf < 1} 

The map G c is defined by the formula G c (x, y) = (£, rj) where £ and 77 are given by equations[T]- 
[2j For better understanding, one should consider c a parameter, and think of the mapping 
G : U ->■ R 2 defined by G(£, T], c) = G c (x, y), where U C R 3 : 

U = {{x,y,c) eR 3 : < c < 1, < x 2 + {y - c) 2 < 1} 

Occasionally, there is a technical advantage to including c in the set of variables, for instance, 
when stating joint continuity, differentiability, etc. which includes the parameter. 

Geometrically, G c acts on a point Q(x,y) in an an almost obvious way. The preliminary 
idea is to map it to the point Q(x,y) to the point P(£,rj), where £ and f] are computed 
from equations [Tp] However, this point is subsequently identified with a point Qi(x',y') 
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where numerically x' — £ and y' = rj. Thus, the action of G c on a point Q is described as a 
composition of two maps ("walks"): 

where the two walks may be described as fallows: 

1. We start at Q, and walk towards S and pass it, until we have walked a total distance 
of 1; 

2. We walk from P towards O and pass it, until we reach Q 1: and satisfy the distance 
condition \QO\ = \OQi\. 

In short: 

G C (Q) = Qi 

For every c G (0, 1) the domain of G c is the open disk about S of radius 1, without its center 
(a punctured disk). The reader should note that according to our conventions A and S like 
on the same side of O when c < 1/2 and on the opposite sides when c > 1/2. Although it 
is possible to enlarge the domain with some additional conventions, we will not do so, and 
adhere to the "natural" domain described above. 

The role of the substitution c i— > 1 — c is explained in the following 

Proposition 1. For an arbitrary c G (0, 1), Let Q be in the domain ofG c and le Q\ = G C (Q). 
Then Gi_ c is well defined at Qi and 

Gi-c(Qi) = Q- 

In particular, the mappings G c and Gi_ c are inverses of each other and G C (U C ) = U\- c . In 
addition the mapping G c : U c — )■ U\- c is a diffeomorphism. 

Proof. Let us consider point P(£,r)) and another point, Pi (£,fj) defined as the reflection of 
Q(x,y) through the point O. We claim that 

1. Pi, R and Q 1 are collinear; 

2. |PiQi| = 1. 

Then the equation Gi_ c (Qi) = Q follows from the above claims and the definition of G in 
terms of walks. Indeed, the claims imply that the two walks defining G\- c {Qi) are: 

Qi h> Pi ^ Q. 

Both properties follow from the observation that the quadrilateral QPPiQi is a parallelo- 
gram, Indeed, its two diagonals are bisected by O. In particular, the side P\Q\ is parallel to 
PQ which has length 1. Therefore, both sides have length 1. Point S lies on the side PQ by 
definition. Thus, R lies on the opposite side P\Q\ because O bisects RS by definition. □ 
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This above statement and proof carefully avoid complicated algebraic equations. An 
attempt to prove the above proposition by brute force calculations is likely to fail. For a 
reader wanting to understand an algebraic approach to the above lemma, we included a CAS 
program Appendix C which arithmetically verifies the claims in the above proof. 

The significance of Proposition [T] in the context of Helfenstein's paper: the invariant sets 
of G c and G^ 1 = Gi_ c are identical. It should also be noted that the following two conditions 
are equivalent for a fixed c: 

1. A function y = f(x) satisfies the functional equation |3j 

2. The graph Cb = {{x, f c {x))} is an invariant set of G c . 

The invariance should be understood locally in the neighborhood of B or near x = 0. Local 
invariance of a Jordan arc Cb means that Cb is contained in the domain of G c . there is a 
neighborhood K of B such that C B H K = G C {C B ) H K. 

The uniqueness of the solutions easily implies that Helfenstein's method cannot work. Let 
us denote by f c (x) any solution of the functional equation |3j defined in some neighborhood 
of x = 0. If we know that the solution is unique then y = f c {x) is a solution of the 
functional equation for both c and d — 1 — c. Uniqueness implies f c {x) = f\- c (x) and 
equality a n (c) — a n (l — c) follows for all n. 



4 Local existence and uniqueness 

It turns out that a properly formulated existence and uniqueness theorem eliminates Helfen- 
stein's approach as viable, but it also eliminates the possibility that a continuous 2e-curve 
exists. The key to such a strong result is a consideration of solutions to the functional equa- 
tion [3] which do not satisfy /(0) = 0, but instead /(0) = yo, where y varies in the range 
1 2/o I < min(c, 1 — c). Such a family provides a good coordinate system near B. Figure [3] 
schematically depicts the situation. In its caption, the figure states several claims which will 
be shown to reflect what has actually been proved about the functional equation and its 
solutions. 

This section is mainly expository, as the proof can be extracted from our Inventiones 
article pQ. The technique is covered in the monograph of Hirsch, Pugh and Shub [4J. There- 
fore, we walk the reader through the constructions and provide some motivations leading up 
to the theorem, which we formulate at the end. 

Throughout this section we use the following set 



V=\(y ,c)eR 2 : 0< 



1 

c 

2 



1 . 



We shall consider a family of curves in IR 3 , {r(y ,c)}, locally represented as a graph y 
F(x,y ,c) and passing through point (0,y ), i.e. 

y = F c (x,y ). 
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Figure 3: The invariant family of curves y = F(x,yo,c) near B for the map G = G c . We 
assume that 1/2 < c < 1/2, which implies that R is between and B. Only a half of each 
curve is drawn in the upper halfplane above the line SB. A typical invariant curve C such 
that C = G C (C) is depicted, as it enters the neighborhood of B foliated by curves of the 
family. A sample trajectory = Q, Q', Q", Q'", . . ., where i = 0,1,2,..., under iteration 
of G, is depicted. Unless even and odd subsequences and Q( 2l+1 \ which must converge, 
both converge to the same point B, the invariant curve C must oscillate between two points 
lying on the straight line RS. 

When c is fixed, we use the notation F c (x, yo) instead of F(x, yo, c) and T c (y ) for the curve 
in M 2 given by the equation y = F c (x,yo). This emphasizes c's role as a parameter. 

We allow (yo,c) G V; equivalently, when c is fixed, we require that |y | < min(c, 1 — c). 
Moreover, we require the local invariance condition: the curve T c (y ) is mapped by G c to 
another curve of the family. Which one? It is easy to verify that 

G c (0,y) = (0,-y). 

Thus, we know that the point (0, — yo) lies in the image of T c (c) and thus 

G c (T c (y )) = T c (-y ). 

locally in a neighborhood of x — 0. 

It can be seen that this invariance condition is equivalent to a functional equation that 
the function F(x,y , c) must satisfy at least in a neighborhood of the line x = 0: 

tp ( x \ F(x,y ,c)) 

y^/x 2 + (c- F(x,y ,c)) 2 J ^x 2 + (c- F(x,y ,c)) 2 

Conceptually, both y and c are parameters in F and once we fix them, we also consider the 
function F y0iC : J — > R, J C M, given by: 

F voA x ) = F ( x ,yo,c). 



n 



Also, conceptually, F is a 2-parameter family of ordinary 1-variable functions. We will use 
the resulting three notations for F interchangeably. The domain J of F y0tC shall depend on 
the parameters. We will assume that J is an open interval, symmetric about x — 0, and 
that its length varies continuously: 

J = J v (yo, c) = {x ER : -<p(y , c) < x < ip(y , c)} 

where if : V — > R is a certain positive, continuous function. Thus, equivalently F : — > R 
where J v C R 3 is given by: 



J v = Ux,y ,c)ER 3 : 0< 



< J> \Vo\< min(c, 1 - c), |x| < </?(y , c) 



2 

For fixed c, we will also consider the section of J' ip : 

J<p(c) = {(x,y ) E R 2 : 1 2/o I < min(c, 1 - c), \x\ < <p(y , c)} . 

Let 

^fGc ifol/2, 
c [G- 1 if c < 1/2. 

We recall that G^ 1 = G\_ c . The reason for this definition is that the theorem below is 
applicable either to G c or G^ 1 depending on whether c > 1/2 or c < 1/2. To maintain 
symmetry of these two cases, we formulate the theorem for H c instead of G c . We recall that 
G^ 1 = Gi- C which means that it would be sufficient to consider only the range c > 1/2. 

Theorem 1. There exists a continuous, positive function tp : V — > R such that 

1. The set J^ic) is contained in the domain of H c . 

2. We have H c (J v (c)) C J v (c). 

3. For every (x,y) G J^c) the limit lim^oo H™(x, y) exists and it lies in the set 

V c = {(x,y) : x = 0,\y\ < min(c, 1 - c)}. 

4- The mapping n c : J$(c) — > V c defined by 

n c (x,y) = lim H£(x,y) 

n— >oo 

is real-analytic and it is a fiber bundle. In particular n c (0,y) = (0,y); equivalently 
tt c \V c = id Vc . 

5. The sets T c (y ) defined for all y s.t. \y \ < min(c, 1 — c) by the formula: 
T c (y ) = n-^yo) = Ux,y) E J v (c) : lim H n c (x,y) = (0,y )| 
may also be equivalently defined by the equation 

V = F(x,y ,c) 
where F : — >■ R is a unique, real- analytic function. 
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6. For all y s.t. \y Q \ < min(c, 1 — c) we have 

H c {T c {y ))cT c {-y ). 

In particular, 

H 2 C (T c (y )) C T c (y ) 

i.e. the curve V c (yo) is invariant under the action of FL\. The curve r c (0) is invariant 
under H c . 

7. We have the following decomposition: 

Mc) = |J T c (y ). 

Vo ■ l?/ol<min(c,l-c) 

8. The convergence in the definition ofT c (y ) is exponentially fast. More precisely, there 
exist continuous functions /i, D : V — > K such that < /x < 1 and < D < oo and 
such that for all (x,y) G ^(c) and aW nonnegative integers n: 

\\H2(x,y) - (0, 2/ )|| < D(y ,c)ti(yo,c) n . 

9. Each curve T c (yo) is tangent to the x-axis, i.e. it is normal to the line RS . Alterna- 
tively, for every (yo, c) G V: 

dF(x,y ,c) = 

An outline of the proof is a subject of another section. Let us note that the curve T(0, c) 
and the corresponding equation y = F(x, 0, c) solve Helfenstein's functional equation [3j 



5 Consequences of local existence and uniqueness 

We shall start with perhaps the least understood consequence of Theorem [T] the non- 
existence of non-differentiable, continuous 2e-curves. We will show only local differentiability 
here, near points A and B. 

Proposition 2. Let c G (0,1/2) U (1/2, c). Let K C M 2 be a closed subset contained in U c , 
the domain of G c , such that G C (K) = K . Moreover, let S ^ K . Let if, J7^(c) and tt c be given 
by Theorem^ Then 

EevnK 

Proof. The assumption S ^ K implies that K is a compact subset of U c . Indeed, otherwise 
it would either have S in it, or contain a sequence Q n G K converging to the circular portion 
of the boundary of U c . Invariance implies that G(Q n ) converges to S, so S G K, which 
contradicts our assumptions. 

To prove the main claim, let us consider a point Q G J^{c) fl K. Let us consider the limit 
E = lim^oo H™(Q) which exists by Theorem [TJ The set K is closed and invariant under 
H c , and thus E G K. Hence, Q G k~ x {E). □ 
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Proposition [2] can be applied to a hypothetical 2e-curve C of G c . Since C must intersects 
the line RS at exactly two points A and B then K must be contained in the union of two 
analytic Jordan arcs tt~ 1 ({A, B}). So, C is automatically locally analytic near points A and 
B. 

Proposition [2] may also be used to deduce an interesting property of the real (not- 
hypothetical!) curves depicted schematically in Figure [2] It depicts two G c -invariant curves 
which are reflections of each other and widely oscillate near one of the points A or B. Let 
C be one of these curves. The curve is an immersed copy of K, but it is not a submanifold 
of M 2 and is not closed. Nevertheless, we may apply Proposition [2] to the closure of C. We 
obtain the following result: 

Proposition 3. If C is an invariant set of G c which is a submanifold disjoint from the set 
x = 0, and let C be its closure. Then either C is an analytic curve or contains an interval 
which is a subset of the line x = 0. 

The interpretation of this result is that either the curve can be completed to an analytic 
submanifold of the plane by adding to it the points of the line x = which lie in its closure, 
or it must widely oscillate near the line x = 0. This fact is interesting because there have 
been many attempts at constructing a 2e-invariant curve by starting with a point Q and 
its image Q' = G(Q), and then connecting them somehow by a path. One can obtain an 
invariant curve C by iterating this path. Authors of these attempts often make a claim of 
continuity of their curve as it approaches the line RS. The above proposition shows that 
such attempts must fail to construct a continuous solution, unless the constructed curve is 
analytic near A and B. As we will point out, the Inventiones article [1] proves that if C is 
a 2e-curve which is locally analytic near the line RS then it is globally analytic. The article 
disproves the existence of analytic 2e-curve, and thus locally analytic also, and furthermore 
by above proposition, continuous 2e-curves. 

To prove that local differentiability of C near A and B implies global differentiability 
at all points, a global proof is required. Such a proof is given Theorems 3 and 4 in the 
Inventiones article [lj. We refer the reader to the proofs there. Let us just comment that 
to prove global differentiability of C we need two facts. The first fact that the map Gf 1 is 
defined on and differentiable on the entire hypothetical 2e-curve C. The second fact is that 
for every point Q e C the limit lim^oo G ±n (Q) exists and is either A or B. Then the local 
differentiable structure near A or B may be "transplanted" by a suitable iteration of a local 
diffeomorphism G to any point of C. 

Theorem [l] and Proposition [2] imply that for every c G (0,1), c 7^ 1/2, the function 
f{x) = F(x, 0, c) is the unique continuous solution of Helfenstein's functional equation[3j and 
this solution is automatically analytic. A simple consequence of this uniqueness is that every 
solution is automatically even: f(x) = f(—x). Indeed, if f(x) solves the functional equation 
then so does f(x) = /(— x), by inspection. Hence, in view of uniqueness, f(x) = f(x). 

To gain a little bit more clarity, we make the following definition: 

Definition 2. We will call a pair of continuous Jordan arcs (Ca,Cb) a local equichordal 
configuration with respect to the points R and S iff: 
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1. The intersection of any straight line parallel to RS with Ca or Cb consists of at most 
one point. 

2. The intersection of Ca and Cb with the straight line RS consists of exactly one point, 
denoted A and B respectively. 

3. For every pair of points P and Q such that: 

(a) P is in Ca and Q is in Cb! 

(b ) the points P and Q and one of the points R or S, are collinear; 
the distance \PQ\ = 1. 

4- The arcs Ca and Cb are symmetric with respect to the reflection in the point O, the 
center of the segment RS , and with respect to the axis RS. 

Less formally, if we only look at the chords of C whose ends are close to A and B, points R 
and S apppear equichordal for the curve C. We request that a local equichordal configuration 
have the symmetry properties, which would follow if a 2e-curve existed. However, we do not 
require that a local equichordal configuration be a part of a 2e-curve. 

It is clear to a reader of Helfenstein's paper that he attempts to prove non-existence 
of a local equichordal configuration and deduce from it that a 2e-curve does not exist. 
The twist is that there exists a local equichordal configuration, but still no 2e-curve exists. 
To Helfenstein's credit, in 1956 figures like Figure [2] were uncommon, used primarily as 
counterexamples in topology. The classical example is the curve 

1 

y = sin — 

x 

which, together with the y axis, joined somehow together to make the figure connected, 
form an example of a set that is not locally connected. Later on, the science of chaos made 
a discovery that curves like this are common in studying differential equations describing 
real mechanical systems. A search on the terms homoclinic connection and heteroclinic 
connection yield many references to such systems. It should be noted that first examples of 
this kind were known already to Poincare and thus available to Helfenstein. 

The reader of Helfenstein's paper can easily verify that Helfenstein tries to prove the 
non-existence of f(x) described by Theorem [TJ Since Theorem [l] shows that the localized 
version of the equichordal point problem has a solution for all c in 0, 1, c ^ 1/2, this means 
that any local non-existence argument focused on a small neighborhood of the line RS must 
fail. This is a deeper reason why the problem had remained open for 80 years until our 
article pQ. 

Corollary 1. For every c in the range < c < 1 Helfenstein's functional equation has a 
unique solution f(x) defined in a neighborhood of x = 0, and which is infinitely differentiable 
and for all n. Moreover, the solution for c and 1 — c is the same, which implies that for all 
c: 

a n (c) = a n (l - c). 
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This corollary is the consequence of the symmetry with respect to the straight line passing 
through O and perpendicular to the line AB. This symmetry exchanges the role of c and 
1 — c. Thus, if a function f(x) solves the problem for c then it also solves it for c replaced 
with 1 — c. Helfenstein was clearly aware of this corollary, but made an incorrect assumption 
about the existence and uniqueness. 

6 Equichordal point problem solved, after all 

Although a n (c) = a n (l — c) for all n, the absence of 2e-curves is proved by different means in 
the Inventiones article pQ. Thus, the line of reasoning used by Helfenstein proved to be one 
of the numerous traps that await anyone studing the problem. The Helfenstein's invariance 
condition appears to have no significance in the Inventiones solution. 

As we have shown, the function f(x) with all the symmetry properties, which is also 
a solution to the Helfenstein's functional equation, does exist locally in a neighborhood of 
x = 0. It is the locality of the argument that is the main error in Helfenstein's paper. Only 
a global argument, which takes into account the behavior of f(x) to the point of breakdown 
of its properties as a single-valued function, can eliminate the possibility of a solution to the 
equichordal point problem. 

Helfenstein's argument is purely local, i.e. it does not refer to the behavior of the function 
outside a small neighborhood of zero. It is clear that the entire oval cannot be represented 
as a graph of a single- valued function f(x). For instance, the standard unit circle is often 
represented as the graph of the equation y = ±\/l — x 2 , where \x\ < 1. However, the 
right-hand side has two possible values. Obviously, any convex oval can be represented by a 
two- valued function with two "branch points", when the two values coalesce to form a closed 
curve. 

Thus, it is clear that the breakdown of the representation of C as a graph of an equation 
y = f(x) must occur at some point, something that Helfenstein does not discuss. Other 
papers deal with this issue. Wirsing in his 1958 article [7], and Shafke and Volkmer in their 
1992 article [3], represent the curve in polar coordinates and use the equation r = g{9) 
to represent the curve. This equation does not suffer from the limit on the range, and is 
capable of capturing the solution to the Helfenstein's equation near A and B simultaneously. 
Moreover, g{9) = 1/2 represents the unit circle, which naturally plays a special role in the 
asymptotic considerations as c — > 1/2. The Wirsing and Shafke and Volkmer research reveals 
the nature of the obstruction to the existence of an oval with two equichordal points: if a 
solution is well-behaved near A, it must lose continuity near B and vice versa. Hence, there 
is no globally defined solution in polar coordinates, either. 
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7 An outline of the proof of Theorem [T] 

The first step of the proof it to interpret the solution of the functional equation [6] together 
with its side condition: 

z \ ~(c- F(x,y ,c)) 

x,-yo,c\ = . = = + (1 - F(x,y Q ,c)), 



'X 2 



(c- F(x,y ,c)) 2 ' J y / x 2 + (c- F(x,y ,c)) 

Vo = F(x, y , c 



as a question about invariant manifolds. We define a map 

G ■ U ->• M 3 , where [/ C M 3 has 

already been defined before: 

U = {(x,y,c) 6l 3 : < c < 1, < x 2 + (y - c) 2 < 1} . 

defined by the formula Q(x, y, c) = (£, r], c) where 

c x 

y/x 2 + (c - y) 2 

v = r ic ~ y) +(i-y). 

a/x 2 + [c-y] 2 

Also see equations [l]|2} It is also true that ^(x, y, c) = (G c (x, y),c), using our prior notation, 
i.e. it simply extends G c to three dimensions by adding a trivial action on the parameter c. 

It can be seen that the surface W given by the equation y = F(x,yo, c) is locally invariant 
under the mapping Q iff F(x,yo,c) solves its functional equation, i.e. Q{W PI U) C VF. 
Because |y| < min(c, 1 — c), we do not have to worry about issues such as non-differentiability 
of Q(x,y,c) or Q~ l (x,y,c) when x 2 + (c — y) 2 = 0, i.e. at the point (0, c, c); it is outside of 
V. Moreover, every point of V is a periodic point of of period 2, because it is easy to see 
that: 

Q(0,y,c) = (0,-y,c). 

Thus, the point (0, 0, c) (corresponding to the point A in Helfenstein's paper) is a fixed point 

ofg. 

In particular, V is an invariant manifold of dimension 2. It is a union of two open 
segments of a straight line. The next step in the proof is to notice that the manifold V is 
normally hyperbolic in the sense of Hirsch, Pugh and Shub [4j. Without repeating lengthy 
definitions, we will explain what this means. We start with linearizing Q at all points of V: 

G(u,y + v,c + w) = G(0,y,c) + DQ(0,y,c) ■ (u,v,w) + 0(\\(u, v, w)\\ 2 ). 

where DQ(0,y,c) is the (Frechet) derivative of Q, i.e. a 2 x 2 matrix. With a little bit of 
work, or using Maxima program presented in Appendix B, we find: 



DG{0,y,c) 



-^--1 

y—c 

-10 

1 
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It happens that DQ(0,y,c) is diagonal, and thus it has three eigenvectors: ei = (1,0,0), 
e 2 = (0, 1, 0) and e 3 = (0, 0, 1), with eigenvalues: 



Ai 

A 2 
A 3 



c-y 



respectively. 

What is important is that the tangent space to V, spanned by the set {e 2 , e 3 }, is invariant 
and |A 2 | = |A 3 | = 1. This means that under iteration of Q the manifold V does not expand 
or contract. It is also important that in the normal direction to V we have linear contraction 
or expansion. This can be verified by performing two iterations of Q, starting from the point 
(0,?/, c). In two iterations, the direction ei scales by the product of the first eigenvalues for 
DQ(0,y,c) and DQ(0, — y,c), i.e. the multiplier 

\c-y J \c + y J c 2 -y 2 

In particular, for \y\ < A/min(c, 1 — c) the multiplier /i is positive. If < c < 1/2, the 
numerator is larger than the denominator and [L < 1. Ifl/2<c<l then the numerator is 
smaller than the denominator, and \l > 1. Of course, when c = 1/2 then /i = 1. 
Let us consider the partition V = V + U V~ , where 

V + =Vn{(x,y,c) : of}, 
V- =Vn{(x,y,c) : c< |} . 

The manifolds V+ and V- are both invariant and satisfy the definition of normal hyperbolicity, 
i.e. under the linearization the normal direction (which in our case is ei) is scaled by a 
multiplier \i ^ 1. The multiplier fi may depend on the point of V± (and it does!), but it may 
not approach 1 along the trajectory of a point, i.e. along the set 

{(0, y, c), (0(0, V, c), G(G(0, y, c)), 0(0(0(0, y, c))), . . .}. 

Since in our case this set consists of two points: 

{(0,3/, c), (0,-1/, c))}. 

Thus, in our case of normal hyperbolicity is easy to verify. 

The theory of normally hyperbolic invariant manifolds immediately implies Theorem [T] 
In passing, we verified all technical assumptions necessary to apply it. Stating the detailed 
results would involve a large amount of notations and definitions, so we shall not repeat it 
here, and refer an interested reader to the classical presentation j3]. 
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We note that our setup satisfied the definition of "immediate absolute r-normal hype- 
bolicity" and thus it satisfies the strongest assumptions used in the classical monograph of 
Hirsch, Pugh and Shub [3]. Thus, the strongest conclusions also hold. It takes very superfi- 
cial understanding to conclude the existence of a continuous F(x, c), infinitely differentiable 
over x, with derivatives continuous in c. This suffices to prove a n (c) = a n (l — c). With 
better understanding, one can show joint analyticity in x and c. 

It should also be noted that our presentation reflects the modern view point, but Wirsing 
[7] proved a theorem in Abschnitt 4 similar to Theorem [Tj He covers only the analytic case, 
but the proof is very concise and to the point, and even today it has some appeal despite 
the lack of generality. 

8 Conclusions 

By repeating the calculations in Helfenstein's paper |2j we identified the mistake in the paper. 
The existence of the mistake was asserted by Wirsing [7j, but he did not directly identify the 
place in which it occurred, which led to half a century of confusion in regard to the validity 
of Helfenstein's paper. We resolved this issue to our satisfaction. In addition, we showed 
that Helfenstein's method cannot be a basis of a new, simpler solution of the equichordal 
point problem than the one currently known [T]. 

9 Notes and references 

9.1 Theorem Q] 

Theorem [I] corresponds to Theorem 2 of the Inventiones article p]. There is a difference of 
notation, because in these notes we adopted Helfenstein's notation. The symmetry properties 
are introduced in the text of the Inventiones article in the parts leading up to Theorem 2. 

It should be noted that Theorem [T] is in the "easy" part of the article, and therefore any 
determined reader is capable of understanding it with minimal effort, given some familiarity 
with invariant manifold theory. Today, invariant manifold theory is quite well understood, 
with many excellent expositions. We still prefer the work of Hirsch, Pugh and Shub [3]. 

The analyticity of F(x, c) follows from an argument of [I] used to prove Theorem 7 there. 
However, the argument is not easy to separate from a more complex situation it is addressing. 
Finally, the analytic version of Theorem [T] follows from Abschnitt 4 of Wirsing's 1958 article 
[7J. Wirsing was not the first one to invent the method of proof. It goes back to the works 
of Hadamard and Perron, and this fact is well known today. The method is used by Hirsch, 
Pugh and Shub and the references to Hadamard and Perron can be found there. 

9.2 Comments by B. Griinbaum 

The following comments can be found in the notes by B. Griinbaum dated 2010 included as 
part of not yet published revision of the book by V. Klee [5]: 
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But more unexpected has been the reaction to the work of Helfenstein [Hel56]. 
His rather infelicitously formulated claim is: "We shall show in this paper the 
non-existence of real and, in one special point, at least six times differentiable 
2e-curves." What his argumentation shows (assuming there are no errors, and 
nobody pointed out any specific errors in the paper) is that a contradiction is 
reached from the assumption of existence of a 2e-curve together with the assump- 
tion that the curve assumed to exist is six times differentiable at a specific point. 
This would seem a reasonable result, and the non-existence would be established 
provided the differentiability assumption could be proved for all 2e-curves as- 
sumed to exist. Hence one would think (and this was explicitly stated by Suss 
[Siis]) that the non-existence of u u 2e-curves would follow from Helfenstein's re- 
sult together with the following result of Wirsing [Wir58]: If a 2e-curve exists, it 
must be analytic, that is, infinitely differentiable at all points. Indeed, to belabor 
the completely obvious, if a curve is infinitely differentiable at all points, then it 
is six times differentiable at a specific point — and the contradiction reached by 
Helfenstein proves the non-existence. It is mystery to me why Wirsing thought 
that Helfenstein's work must be in error since "it is contradicted by the present 
investigation" ("durch die vorliegende Untersuchung widerlegt wird"). But an 
even greater mystery is why Klee, in all his discussions of equichordality, ac- 
cepted Wirsing's statement as valid, and Helfenstein's result as invalid. Wirsing 
committed a logical error, and Klee and others uncritically accepted it a valid. 
The reviews of Helfenstein's paper in the Math. Reviews and the Zentralblatt 
fail to claim any errors in it. Not all details of Helfenstein's proof are given in 
[Hel56] - but no error has been found either. 

Of course, [Hel56] refers to the earlier cited article [2J. The Wirsing's work cited as [Wir58] 
is [7]. The last citation in the above comment, [Suss], is [5]. 

Clearly, this inaccurate account of the state of the equichordal point problem still lingers 
in the public domain. Griinbaum writes in the introduction ([5], p. i): 

About mid-May 2010 it occurred to me that it might be appropriate to have 
the fifty-years old collection made available to participants at the "100 Years in 
Seattle" conference. 

In contrast with Griinbaum, we are certain that Wirsing rejected Helfenstein's argument 
for the right reasons: it contradicted Wirsing's own research. Wirsing refers to Helfenstein's 
work in two places. The main point is made in his Abschnitt 4, in which Wirsing proves a 
variant of our Theorem [TJ and at the end he writes: 

Es bleibt die Frage, ob fiir irgendwelche c und X die analytische Fortsetzung von 
€\ and €2 zu geschlossen Kurven fiihren kann. 

Ubrigens steht die Tatsache, dafi die Doppelspeichenbedingung sich wenigstens 
lokal (in Umgebungen von T\ and T2) durch regular analytische Kurven befriedi- 
gen lafit, im Gegensatz zu der Arbeit von Helfenstein [5], der durch 6-malige 
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Differentiation der Funktionalgleichung bei T\, T 2 zu einem Widerspruch gelangt 
und daraus auf die Nichtenxistenz einer 6-mal differenzierbaren D-Kurve schliefit. 
Der Fehler diirfte in den unveroffentlichten Rechnungen liegen. 



10 Asymptotic analysis beyond all orders 

Shafke and Volkmer [3] conducted a deep asymptotic analysis of the equichordal point prob- 
lem, which implies that a 2e-curve may exist only for a finite set of values of c. In spirit, 
the paper continues the researches of Wirsing. Wirsing proved that if c — > 1/2 then the 
hypothetical 2e-curve must be extremely close to the circle of radius 1/2 centered at O, 
closer than C a |c — l/2| a where a is an arbitrarily large power, and C a is a certain constant, 
depending on a. Thus the name "asymptotic analysis beyond all orders" which is sometimes 
used in reference to this kind of result. 

Shafke and Volkmer quantified this statement even further, expressing the leading asymp- 
totics of the difference between the 2e-curve and the circle in terms of the exponential. In 
fact, the paper extends the representation of the arcs Ca and Cb in the discussion following 
Theorem [I] to the full angle 9 in polar coordinates except one point (either zero or it) near 
which the arcs start rapidly oscillating, breaking continuity! 
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A Maxima code showing Helfenstein's mistake 



/* _*_ Mode: Maxima -*- 

* Author: Marek Rychlik, Date: January 20, 2012 

* This program finds the coefficients a[n] of the Helfenstein's 

* paper by expanding f (x) in a Taylor series about x=0 

* and comparing coefficients at like powers. 

* Additionally, it factors a[n], and finds b[n] obtained 

* from a[n] by substitution c = 1/2 + sqrt(z)/2. 
*/ 

assume (c>0) ; 

/* Define Taylor expansion of f (x) with generic coefficients a[k]; 

the meaning of a[k] is consistent with Helfenstein */ 

def taylor (f (x) , sum (a [k] *x~k,k, 2, inf ) ) ; 

/* The denominator */ 

d(x,y) :=sqrt(x~2+(c-y)~2) ; 

ksi(x,y) :=x/d(x,y)-x; 

eta(x,y) : =- (c-y)/d(x,y)+(l-y) ; 

/* The functional equation */ 

eqn: f (ksi (x,f (x) ) )=eta(x, f (x) ) ; 

/* The maximum derivative order we consider */ 

nlimit : 10 ; 

/* Expand functional equation in a Taylor series */ 

tay_eqn : t ay lor ( eqn, x,0, nlimit) ; 

/* Extract coefficients at relevant powers */ 

sys_eqn:makelist (coef f (tay_eqn,x,k) ,k,0 , nlimit) ; 

/* Solve for the coefficients a[k] */ 

soln: solve (sys_eqn,makelist (a [k] , k, 2, nlimit ) ) ; 

/* Assign coefficients */ 

for n:2 thru nlimit do ( 

a [n] : ev (a [n] , soln) , 

a_f actored[n] : factor (a [n] ) 

); 

/* Translate coefficients to TeX */ 
for n:2 thru nlimit step 2 do ( 

tex('a[n] = a[n]) , 

tex('a[n] = a_f actored[n] ) 

); 

/* Find the coefficients b[k] */ 

for n:2 thru nlimit step 2 do b[n] :ratsimp(ev(a[n] ,c=l/2+sqrt(z)/2)) ; 
/* Translate coefficients b[k] to TeX */ 
for n:2 thru nlimit step 2 do ( 
tex('b[n] = b[n]) 

); 
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B Maxima code finding Frechet derivative of G 



/if _*_ Mode: Maxima -*- 

* Author: Marek Rychlik, Date: January 20, 2012 

* This program finds the derivative DG(x,y,c) 

* and simplified expression for DG(0,y,c). 
*/ 

kill(all); assume(c>0); assume(y<c); 
/* The denominator */ 
d(x,y,c) :=sqrt(x~2+(c-y)~2) ; 
ksi(x,y,c) :=x/d(x,y,c)-x; 
eta(x,y,c) : =- (c-y)/d(x,y, c)+(l-y) ; 
/* The mapping G */ 
G(x,y,c) :=[ksi(x,y,c) ,eta(x,y,c) ,c] ; 
/* The Frechet derivative of G */ 

A : matrix (diff (G(x,y , c) ,x) , dif f (G(x,y , c) ,y) , dif f (G(x,y , c) , c) ) ; 
/* Typeset the derivative as TeX */ 
tex(A) ; 

/* Compute the Frechet derivative at a point of V */ 
B:ev(A, [x=0] ) ; 

/* Clean up by expanding into partial fractions */ 
B:partf rac(B,y) ; 

/* TeX form of the derivative at a point of V */ 
tex(B) ; 



C Maxima code to verify formula for inverse of 

/* _*_ Mode: Maxima -*- 

* Author: Marek Rychlik, Date: February 1, 2012 

* This program provides an arithmetical proof that 

* G_c and G_{l-c} are inverses of each other. 
*/ 

/* The denominator */ 
kill(all) ; 

assume (c>0 , c<l , y<c , y-c+l<0 , y+c>0 , y+c- 1>0) ; 

/* The denominator */ 

d(x,y,c) :=sqrt(x~2+(c-y)~2) ; 

u(x,y,c) :=x-x/d(x,y,c) ; 

v(x,y,c) :=y+(c-y)/d(x,y,c) ; 

ksi(x,y,c) :=-u(x,y,c) ; 

eta(x,y,c) :=l-v(x,y,c) ; 

/* The map */ 

G(x,y,c) :=[ksi(x,y,c) ,eta(x,y,c)] ; 
/* The generic point near B */ 
Q: [x,y] ; 
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/* Point P is collinear with S and Q, such that dist(P,Q)=l */ 
/* Note: P is expressed in xy-coordinates */ 
P: [u(x,y,c) ,v(x,y,c)] ; 

/* Reflection of Q; also expressed in xy-coordinates */ 

Pl:[-x, 1-y]; 

/* Ql = G(Q) */ 

Ql: [ksi(x,y,c) ,eta(x,y,c)] ; 

/* Checks that dist(Ql,Pl)=l */ 

/* Must print "true"! */ 

is(ratsimp(apply("+", (Ql-Pl) ~2) )=1) ; 

/* Checks that P1,Q1 and R are collinear (ksi-eta coordinate system) */ 
R: [0,1-c] ; 

/* Must print "true"! */ 

is (rats imp (determinant (matrix (Pl-R, Ql-R) ) )=0) ; 
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